Intrinsic spin-orbit interactions in fiat and curved graphene nanoribbons 
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Recent theoretical and experimental works on carbon nanotubes and graphene samples have 
revealed that spin-orbit interactions, though customarily ignored in carbon-based materials, are 
more important and complex than it was thought. We study the intrinsic spin-orbit coupling effects 
on graphene nanoribbons, both flat and bent. Calculations are performed within the tight-binding 
model with the inclusion of a four-orbital basis set. Thereby the full symmetry of the honeycomb 
lattice and the hybridization of a and it bands are considered. In addition to the zero-energy tt- 
edge states, cr-derived edge states are found for the three investigated ribbon geometries. The a 
states are also spin-filtered and localized at the boundaries of the ribbons. The calculated spin-orbit 
splittings are larger for the a- than for the 7r-derived edge states. Due to this enhancement, spin- 
orbit splittings of the er-states reach values of the order of a few Kelvin. These spin-filtered edge 
states are robust under a — n hybridization and curvature effects. 

PACS numbers: 73.20.-r, 73.22.-f, 73.23.-b 
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I. INTRODUCTION 

The influence that geometry and size have on the 
electronic properties of carbon based materials has 
been already proved. Two dimensional graphene 
sheets, graphene stripes known as graphene nanoribbons 
(GNRs), or carbon nanotubes (CNTs) formed by rolling 
a graphene sheet unto a cylinder, due its unique elec- 
tronic properties mainly to their geometries. Chirality 
and diameter determine the properties of CNTs, while 
the GNRs physics depends on its width and the shape of 
its edgesA. This strong dependence of the electronic prop- 
erties on the geometry arises from the peculiar structure 
of graphene with carbon atoms ordered in a honeycomb 
lattice, divided into two nonequivalent triangular sub- 
lattices, and confers these materials promising potential 
applications in electronics and spintronics. The low en- 
ergy physics of carbon-based materials is governed by the 
states close to the Fermi energy in graphene. They cor- 
respond to the 7r-states near the nonequivalent K and K' 
points at opposite corners of the hexagonal Brillouin zone 
(BZ), the Dirac points, at which the valence and conduc- 
tion bands touch and present a conical energy spectrum. 

Due to the low atomic number of carbon, spin-orbit 
interaction (SOI) was expected to be small. Neverthe- 
less, several theoretical studies were devoted to investi- 
gate SOI effects in graphite 2 and more recently in CNTs 
and graphene 3 - - — . In a k.p scheme and using a pertur- 
bation approach, it was shown that the SOI opens a gap 
in the energy bands crossing at the Fermi level and lifts 
the spin degeneracy for higher energy subbands 3 . In a 
previous work 6 , within a four-orbital tight-binding (TB) 
framework, we have demonstrated the chiral dependence 
of the intrinsic SOI effects, which can not be described 
by a continuum model. Spin-splitting is only possible for 
chiral CNTs due to the lack of inversion symmetry, while 
the spin-degeneracy remains for achiral CNTs. Further- 
more, the recently measured asymmetric splittings of va- 



lence and conduction bands in a CNT quantum dot£ are 
also explained within the same modeli^r— pointing out 
that curvature enhances the SOI strength. In fact, recent 
experimental measurements of the energy shifts caused 
by the SOI, have proved that the SOI effects in CNT 
quantum dots and graphene samples are not as small as 
predicted&i 3 -. Besides, since SOI is the key ingredient for 
the Quantum Spin Hall state, the new phase of matter 
proposed by Kane and MeleiA in graphene, the previously 
considered almost negligible SOI interaction in graphene- 
based nanostructures has risen to be a main topic in the 
condensed matter fielcU£~— . 

The electronic properties of GNRs have been inten- 
sively investigated due to the remarkable behavior of the 
zero-modes edge states 2 ^—. Most of the works have 
been based either on the k.p or the effective 7r-orbital 
TB approaches, the simplest models that capture the 
physics of graphene, although some first-principles and 
LCA02£r— calculations have also been reported. Fur- 
thermore, the effect of the SOI in graphene sheets and 
GNRs with different boundaries is currently attracting 
great interest and, in spite of the large amount of pub- 
lished works, many aspects are still hidden. Some nu- 
merical calculations based on multiorbital models 3 ^, and 
first principles— ^ can be found in the literature but, to 
our knowledge, no investigation of the SOI on flat and 
bent GNRs based on a sp 3 tight-binding model has so 
far been performed. 

We address here the study of the intrinsic spin-orbit 
coupling effects on the electronic properties of GNRs 
with edges of different shapes: zigzag, Klein's bearded 
and armchair edges. We consider flat and bended rib- 
bons to investigate the interplay between the SOI and 
the curved geometry and focus on the topological edge 
states, not only the well known 7r-derived states, but also 
those arising from cr-orbitals. We investigate how curva- 
ture influences the interplay of confinement and SOI in 
determining the topological character of the edge states. 
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An explicit description of edge states requires a model 
that gives the energy bands through the entire BZ. There- 
fore, calculations are performed by exactly solving an 
empirical tight-binding Hamiltonian with a sp 3 basis set. 
The TB model handles equally all points in the BZ, which 
is not possible with the continuum model and is especially 
appropriated with regard to the treatment of the SOI. Be- 
sides, it takes into account the full symmetry of the hon- 
eycomb lattice, hence the threefold rotational symmetry 
and the trigonal warping are included. Nor external po- 
tentials or doping are taken into account, although the 
effect of the a-derived edge states on the behavior of the 
spin Hall conductivity with the chemical potential would 
be discussed. 

In addition to the 7r-edge states we found that a- 
derived edge states are also present for all the boundary 
geometries of the GNRs studied. Further, as expected 
its energy is independent of GNR width and for a large 
interval of k values of the one dimensional BZ the a edge 
sates are located in a gap. Hence, they are likewise per- 
fectly localized along the edges of the ribbon. Our main 
finding is that, when the SOI term is included in the 
Hamiltonian, independent of the GNR termination, the 
fourfold degeneracy of edge states, both of it and a char- 
acter, is lifted except at the time reversal invariant sym- 
metry points r and K of the one dimensional BZ. They 
split into the spin filtered edge-states, which are robust 
upon a — it bands hybridization and, as expected for a 
topological state, are unaffected by curvature. 

The paper is organized as follows: in Section II the 
details of the model are given. In Section III the results 
obtained for the different type of GNRs are described, 
and curvature effects are analysed by comparison of flat 
and curved geometries. Finally, we conclude in Section 
IV summarizing our results and drawing some conclu- 
sions. 



II. MODEL AND METHOD 

A ribbon consists of a stripe of graphene of infinite 
length and finite width. There are two prototypical 
shapes of the GNR edges, zigzag and armchair, they have 
a 9 = 30 angle of difference in the cutting direction and 
the termination of a generic ribbon has a combination 
of the two types. We consider ribbons of three different 
kinds of terminations: zigzag, bearded and armchair, and 
define the width of the GNRs as n, where n stands for 
the number of zigzag lines for the zigzag and bearded rib- 
bons and for the number of dimer lines for the armchair 
ribbon 21 . The corresponding unit cell contains N = 2n 
carbon atoms. Examples of the three types of ribbons are 
displayed in Figure Q] We do not consider here any re- 
construction or relaxation of the edge o 34 ' 35 . Nevertheless, 
chemical and structural modifications of graphene edges 
are promising routes to design GNR based device a 20 ' 25 ' 27 . 

Curvature effects are introduced by bending the ribbon 
along its width. The bending is realized without stretch- 




FIG. 1: Schematic lattice structure of graphene nanoribbons 
with: a)zigzag, b) bearded, and c) armchair edges. 



ing, varying the atom coordinates to form an open cylin- 
der. Different curvatures are obtained by changing the 
angles and diameters of the cylindrical configurations. 

The electronic properties of the graphene nanoribbons 
are calculated from the one-electron Hamiltonian given 
by, 



H — Hq + Hso, 
Hq is the spin-independent noninteracting Hamiltonian, 



*o = ^ + V(r), 

and and Hso is the microscopic spin-orbit interaction 
term 
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where V is the full crystal potential, p the electron mo- 
mentum and S represents the spin operator. 

We calculate the band structure of the graphene 
ribbons by the Slater-Koster— empirical tight-binding 
(ETB) Hamiltonian. A four-orbital, 2s, 2p x , 2p y and 
2p z , basis set is considered in order to include the conven- 
tional on-site approach for the intrinsic SO interaction. 
Within the TB approximation the Hq term is written as, 
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where e a represents the atomic energy of the orbital a, < 
ij > stands for atomic sites of the honeycomb lattice and 
c™j~ and cf s are the creation and annihilation operators 
of one electron at site i, orbital a and spin s, respectively. 

Since the major contribution of the crystal potential, 
VV, is near the atomic nuclei, the intrinsic SOI can be 
accurately approximated by a local atomic term of the 
form: 
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FIG. 2: (Color on line) Energy splitting for the conduction 
(circles) and valence (diamonds) bands of the (61,0) zigzag 
CNT as a function of the value of the spin-orbit coupling 
constant A. The solid(black), dashed(red) lines are a guide to 
the eyes. 

where spherical symmetry of the atomic potential is as- 
sumed, ri is the radial coordinate with origin at the 
atom, Vi the spherical symmetric potential about the 
same atom and L the orbital angular momentum. The 
atomic SO coupling constant A depends on the orbital 
angular momentum L and thus on the atomic orbital. 

For the sp 3 model Hamiltonian, SO coupling occurs 
only among p orbitals and neglecting nearest neighbor 
SO terms, the intrinsic SOI is described as an on-site 
interaction among the p-orbitals. In contrast, in the k.p 
and effective 7r-band models, the atomic SOI is treated by 
perturbation theory as a second- neighbor spin dependent 
hopping terrr&i 

Hso adds diagonal and off-diagonal spin-dependent 
matrix elements to the 8N x 8N matrix Hamiltonian, 
where N is the number of carbon atoms of the unit cell 
and 8 comes from the four-orbital per spin TB basis set. 
Using the raising and lowering angular momentum oper- 
ators, L + = L x + iL y and L_ = L x — iL y , respectively 
and the Pauli spin matrices, the complete Hamiltonian, 
H, in the 2x2 block spinor structure is given by, 

Tj f Hq + XL Z AL_ \ ( , 

H -{ XL + H Q -XL Z )- W 

where the atomic-like spin-orbit term Hso has been 
added to SK-ETB spin-independent Hq Hamiltonian. 
The diagonal terms act as an effective Zeeman field pro- 
ducing gaps of opposite signs at the K and K' points of 
the BZ. 

In the sp 3 TB model the non-vanishing matrix ele- 
ments of Hso couple 2p Xl 2p y and 2p z orbitals. These 
terms have been widely discussed previously, see for ex- 
ample Ref , 3 LS : 8 . Hso induces a — it hybridization and 
does not break the time reversal symmetry, therefore spin 
degeneracy can not be removed on systems with inversion 
symmetry. 

The electronic properties of the GNRs are obtained 
by the exact diagonalization of the total Hamiltonian 



H. The Tomanek-Louie parametrization for graphite 39 , 
previously used to calculate SOI effects on CNTs, has 
been also employed in the present calculations. The 
TB parameters, taken from Ref22 are : e s = -7.3eV, 
e p = O.OeU, sso= -4.30eV, spa = 4.98eU, ppa= 6.38eU, 
ppir~-2.66eV, and specify the energy scale of the model. 

The spin-orbit coupling strength of graphene is un- 
known and its exact value is under debate. Due to the 
small atomic number of carbon, it has been assumed to 
be very small and different values in the order of w meV 
have been give n 7 ' 8 ' 14 . In order to have a quantitative es- 
timate of the SOI induced splittings, we have evaluated 
A from the SO induced gaps experimentally measured in 
a CNT quantum doA In the experiment a zigzag small 
gap CNT of about 5 nm of diameter is probed using tun- 
nelling spectroscopy. The anisotropic gaps measured at 
zero magnetic field are of « 0.37 meV for electrons and of 
« 0.22 meV for holes. The (61,0) CNT with a diameter 
of about 4.8 nm is a zigzag semiconducting tube of the 
n = 3m + 1 family, so the SOI induced splitting is greater 
for the conduction than for the valence band 10 as occurs 
in the experiment. Calculated SOI splittings for the con- 
duction and valence bands are represented in FigJ5]for the 
(61,0) nanotube as a function of A. From the figure it is 
inferred that A w 3 meV is needed in order to to obtain 
SOI splittings similar to those experimentally reported. 
This value is of the same order of magnitude of previous 
estimations: 7,8 . In the figures we present results obtained 
with values of A in the range of 0.2 — 0.4eU, for illustra- 
tive purposes. This range corresponds to 0.03 — 0.06 in 
units of the ppa=6.38 eV parameter. 



III. RESULTS 

The electronic properties of GNRs are derived from the 
band structure of graphene subject to a stripe geometry. 
The combination of the confinement due to the finite-size 
and the presence of boundaries yields the peculiar band 
structure of GNRs. The truncation of inter-atomic bonds 
caused by the borders gives rise to the appearance of 
edge states, which are strongly dependent on the atomic 
termination of the GNR. Moreover, the energy subbands 
associated to the intrinsic band structure of the graphene 
sheet are also dependent on the boundary conditions of 
the GNR. Calculations have been carried out in GNRs 
of different geometries and widths. In order to avoid the 
coupling between edge states localized at each boundary, 
ribbons of more than 50 chains are considered. 



A. Zigzag graphene nanoribbons 

In zigzag graphene nanoribbons (ZGNRs) the atoms 
of each edge belong to the same sublattice and oppo- 
site edge atoms are of different sublattices. In the ab- 
sence of SOI, the band structure of ZGNRs presents 
together with the subbands originating from the two- 
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FIG. 3: Band structure of the zigzag ribbon n=60 obtained 
with SOI A = 0.4. k are in units of tt 
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FIG. 4: Zoom of the band structure of the ribbon N=60 with 
A = 0.4 for the two edge states, left around k — 1 for the zero 
edge states (energies are multiplied by 30) and right in the 
region of k — for the a derived band, k are in units of tt 



dimensional graphene electronic structure, zero-energy 
7r-orbital derived states at the |tt < k < tt interval. 
These states, with wave functions localized at the edges, 
form flat bands and give rise to a peak in the density 
of states at the Fermi energy being crucial for magnetic 
instabilitie o 21 ' 40 . Besides, between 1 and 3eV below the 
Fermi energy, a-orbital derived edge band appears in 
the sp 3 TB calculation. These states, which are missed 
in one-7r band models, disperse along the BZ and their 
wave function amplitude is also fully localized at the edge 
atoms around k = and k = ±7r, where the ZGNRs have 
energy gaps. It is worth mentioning that ribbons of dif- 
ferent widths have been studied, from n = 10, to n = 200 
and edge states, both tt and a, appear at similar energies 
independent of the ribbon width. 

To analyze the result of the SOI in ZGNRs, the band 
structure of the n = 60 fiat ZGNR with A = 0.4eV is 
shown in Figure[3]in the range of — tt < k < tt. The spec- 
trum remains gapless and the degeneracy of edge sates, 
originated from both tt- and a- orbitals, is partially lifted. 
A zoom of the dispersion relation of the tt and a edge 
states around the time-reversal invariant k = ±7r and 
r points, respectively, are presented in Figure 01 The 
SOI shifts down in energy the edge states giving a small 
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FIG. 5: (Color on line)Wave function amplitude of the lo- 
calization versus atomic coordinates accross the ribbon width 
for the two edge states states represented in Fig. [4] a) for 
the state at the Fermi level at k = tt (black stars), k — 0.85-7T 
(red circles) and k = 0.7-7T (green diamonds) b) for the a de- 
rived state at k — 0(black stars), k — 0.15-7T (red circles) and 
k — 0.3i7r (green diamonds). Black solid, red dashed and 
green dott-dashed lines are guides for the eyes. 



dispersion to the tt flat bands. It also lifts partially the 
fourfold degeneracy of both tt and a edge states, which 
become twofold degenerate, except at the k = ±tt and 
k = 0/(k = ±7r) points of the BZ, respectively. Each 
edge state splits into two degenerate Kramers doublets, 
with linear dispersion in a very small fc-region around the 
crossing and forward and backward mover states with 
opposite spin. Since degenerate states are confined at 
different edges of the ZGNR, the resulting states with 
opposite spin are localized in either edge. Thus two in- 
dependent spin-filtered edge states are at each boundary 
of the ribbon. Hence, for a given energy the ZGNR has 
four conducting channels spatially separated, an extreme 
of the ribbon contains a forward mover with a given spin 
S and a backward mover with opposite spin, — S, and 
conversely for the other extreme of the ribbon. Note, 
that the energy splitting of the tt and a edge states dif- 
fers almost by a factor of 30. 

The localized character of these states can be observed 
in Figure [31 where the wave function amplitudes are rep- 
resented against the atomic coordinates across the width 
of the n = 60 ribbon, for three different values of k in the 
BZ. The wave function amplitude of both tt and a states 
is fully localized in the edge atoms in a large k interval 
around the k = ±7r and T points, respectively. The lo- 
calization length of the zero-energy state increases with k 
and eventually, for k close to ±t£, where it merges with 
the bulk bands, becomes extended. Further, the spatial 
localization correlates with the tt-u hybridization induced 
by the SOI. At k = ±7r, where the state is composed of 
purely tt orbitals, its wave function is fully localized in the 
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edge atoms, moving towards the T point, i.e. k = ±0.97r, 
the localization starts to decrease toghether with an ad- 
mixture with a orbitals. Therefore, the presence of the 
SOI and consequently of the ir-a hybridization, slightly 
reduces the localization length of tt edge states. Analo- 
gously, the a edge states show an almost pure a orbital 
character and a strong localization at the borders of the 
ribbon for both T and k = ±7r points. The presence of 
a small tt contribution, for k values different from the 
high symmetry points, gives rise to the increase of the 
localization length. 

Due to spin-orbit coupling S is not longer a good quan- 
tum number and the eigenfunctions are a linear combi- 
nation of spin-up and spin-down states. As stated above, 
each extreme of the ribbon contains two states, forward 
and backward movers, with opposite spin S. The cal- 
culated expectation value of < S > shows that the spin 
orientation of both tt and a edge states is almost per- 
pendicular to the graphene plane at the time-reversal in- 
variant k = 7r and T points. The deviation from the 
perpendicular axis is smaller than 1%. The orientation 
axis, although depends on the magnitude of k, slightly 
change in the k interval in which edge states remain spa- 
tially localized. When a-TT hybridization becomes rele- 
vant and states deslocalize, the in-plane component has 
a finite value and the spin orientation axis of extended 
states forms an angle with the graphene plane. Further- 
more, the expectation value of the orbital angular mo- 
mentum < L > for edge states is almost zero, indicating 
the quenching of the orbital angular momentum. 

Finally, because of the estimated small value of A ss 3 
meV , the corresponding SOI splitting of 7r-states is only 
of the order of 0.17 x meV, which gives tempera- 

tures in the range of 10 -3 K, in agreement with previous 
estimations. Nevertheless, for the new a-derived edge 
states the calculated splitting is much larger, 0.51 meV, 
which results in temperatures of the order w 5 K that, 
although small, is experimentally attainable and thus, 
the topological protected spin-filtered states of graphene 
could be observed. 



B. Klein bearded graphene nanoribbons 

In this termination, as in the case of zigzag edges, all 
atoms of one end belong to one sublattice, and those of 
the opposite end, to the other sublattice. The Klein's 
edge could be formed by bonding an additional C atom 
to every edge site of the zigzag ended ribbon^. Zigzag 
and Klein ribbons differ in the number of cut bonds. 

Figure [6] shows the dispersion relation of a n = 60 
ribbon with this geometry without and with the SOI ef- 
fect. The Klein boundary produces as well edge states 
localized at the ends of the ribbon, with the same be- 
havior of those of the zigzag nature, although the zero- 
energy fiat bands lie at the central region of the BZ 2 b 23 , 
— =g- < k < Besides the flat zero edge tt states, there 
are two new edge states of a character. In the absence of 




FIG. 6: Band structure of the bearded edge N=60 ribbon 
without(with) SOI left(right). k are in units of tt 



SOI, one of the a states, not shown in the figure, lies well 
below the Fermi energy, at around -4 eV, mixing with 
the bulk bands, and showing a weak dispersion along the 
whole BZ. The other a state, lies at zero energy at the F 
point, but presents energy dispersion, lowering its energy 
as |k| increases, lying at —0.23 eV for k = ±7r. The flat 
bands are of pure tt character, and are localized in the 
outermost atoms at either edge of the ribbon. Also the 
wave function amplitude of the the dispersive band of 
mainly a nature is strongly localized at the edge atoms, 
i=3 0.91% at the extreme atoms. The spin-orbit coupling, 
besides to induce the splitting into two Kramer's dou- 
blets of the edge bands, has drastic effects on the energy 
and dispersion of both tt and a bands. While tt states 
acquire a small dispersion, a large energy shift and flat- 
tening of the a bands occurs. These effects originate from 
the ir-a hybridization induced by the SOI and depend on 
the value of the coupling strength. For example, the zero 
energy state was at k = 0.27T pure 7r-like in the absence 
of SOI term, but has around 1.1% of a components when 
the SOI term is considered with the A = 3 meV coupling. 
As shown in Figure |5] bottom, the SOI shifts the edge 
bands, inducing a separation at T, for A = 0.4 eV, of 276 
meV which varies almost linearly with the value of A and 
therefore it reduces to 7.38 meV estimated when A = 3 
meV. 

Moreover, the SOI induced energy splitting of edge 
states is of » 0.38 meV and 0.18 x 10 -3 meV for the lower 
and higher energy bands respectively at k — ±0.27r. The 



difference of the splittings is due to the different orbital 
compositions of the bands. The energy splitting of the 7r 
and a edge states is of the same order than in the zigzag 
ribbons. The dependence of the orbital composition on 
k explains the different splitting magnitudes at different 
regions of the BZ. The extended bulk subbands states, 
associated to the intrinsic band structure of the graphene 
sheet are, in the presence of spin-orbit coupling, still dou- 
bly degenerate due to the combination of time reversal 
and inversion symmetries. 

The average < S > and < L > values, analogous to 
those of zigzag ribbons, point out the counterpropagation 
nature of the edge states and give a spin quantization 
axes almost perpendicular to the sheet axes. 



C. Armchair graphene nanoribbons 

Armchair edges are formed by homogeneous lines of 
dimers where atoms belonging to the two different sub- 
lattices alternate. The electronic properties of these rib- 
bons present a strong dependence on their width, as it 
happens for zigzag CNTs. The armchair ended ribbons 
are metallic when its width W = (3M — l)d, where M is a 
positive integer and d is the C-C atom distance. For other 
values of W the armchair ribbons (AGNRs) are semicon- 
ductors with direct band gaps at T which are inversely 
proportional to its width^S^. They follow a family be- 
havior as it is found for zigzag CNT a 10 ' 30 , where metallic 
behavior is obtained in a sequence of period 3. However, 
zigzag CNTs are metallic for n — 3M due to the different 
periodic boundary conditions imposed by the cylindrical 
geometry. Ribbons, with both metallic and semiconduc- 
tor behavior, of different widths are calculated in order 
to clarify how the SOI affects both electronic structures. 

Metallic AGNRs 

Armchair ribbons of width W = (2M — l)d present 
a band structure with two states crossing at the Fermi 
energy at the T point of the BZ. Figure [JJ shows the band 
structure of the n = 62 ribbon calculated with SOI and 
A = 0.4. In the absence of spin-orbit coupling, the cross- 
ing states are only spin degenerate, have pure 7r-orbital 
character, show linear dispersion and their wave func- 
tions are extended throughout the whole ribbon width. 
The inclusion of the SOI term in the Hamiltonian opens 
a small gap at T, for the ribbon n — 62 is of s» 1CP 4 
meV for A = 3 meV, which increases as the width of 
the ribbon decreases^. Further, it changes the orbital 
composition of the linear states allowing for a small pro- 
portion about 5%, of a orbitals and the pure spin nature 
is not conserved. Nevertheless, the two linear dispersing 
states remain twofold degenerate and keep its extended 
character. These states have been found to persist in 
the presence of intrinsic spin-orbit interactions as spin 
filtered states localized on the ribbon edges within the 
low-energy Dirac model^, In the present model, in order 
to obtain localization of these states, we need unphysical 
high values of A of the order of 2.8 eV. 
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FIG. 7: Band structure of the armchair N=62 metallic ribbon 
(top) and N=60 semiconductor ribbon (bottom) with SOI 
A = 0.4. k are in units of n 



Edge states composed mainly of a orbitals and four- 
fold degeneracy appear below the Fermi level. They lie 
at -1.27 eV at the k = ±7r points of the BZ, are disper- 
sive and fully localized at the edge atoms of the ribbon. 
At k w ±0.37r merge with the bulk bands, but keep its 
fourfold degeneracy and localized character. The wave 
function amplitude is confined to the outermost atoms, 
although, at the zone center, localization occurs at both 
ends of the ribbon. These er-orbital derived states are 
truly surface states appearing at the same energy inde- 
pendently of the ribbon width, as we have confirmed cal- 
culating AGNRs of up to n = 122 and n = 242. 

The a derived edge states, upon SOI inclusion, split 
into two doublets, except at the time-reversal invariant 
k = ±7r points, where they remain fourfold degenerate. 
The value of the splitting varies with k reaching a max- 
imun value of « 0.36 meV at k = ±0.57r for A = 3 meV. 
These states keep the localization at the boundary atoms, 
although their orbital composition presents now a small 
proportion of tt orbitals. The so-called bulk subbands of 
the ribbon mantain always their double degeneracy. 

Semiconductor armchair GNRs. 

The bottom of Fig. [7J shows the band structure corre- 
sponding to the semiconductor ribbon N=60, with SOI 
A = 0.4. The energy gap at T is of 0.156 eV almost equal 
to the value obtained without SOI. The a-orbital derived 
edge states and bulk related subbands present the same 
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features as in the metallic ribbons. Even the SOI energy 
splitting of the Krammers doublets are of the same order 
than in the previous case. 



D. Curvature effects 

Next we analize the effect of curvature. As stated 
above, in the curved geometry bond stretching is not 
allowed along the ribbon. Even the distances between 
the carbon atoms at the edges are not changed in any of 
the different terminations studied. The lattice structure 
of the ribbon is isotropically bended along the width, 
x-direction, while the bonding distances, and therefore 
the hoppings, between atoms in the y-direction are not 
modified. 

Curvature induces hybridization between tt and a or- 
bitals. In fact several works have shown that curvature 
enhances the SOI effects in CNTs&±&±i'2I The curvature 
induced a — it hybridization is related with the amplifica- 
tion of the SOI effects for small diameter tubes. However, 
curvature effects turn out to be weaker for ribbons than 
for CNTa2i-i2. The cylindrical shape of the tubes impose 
periodic boundary conditions that confer symmetries to 
the wave functions which do not hold in the open ge- 
ometries. As for the flat geometry, a critical value of the 
width -slightly greater- is needed in order to avoid the 
coupling of edge states in the bent ribbons. Above the 
critical value of the width of the ribbon, the behavior of 
the edge states found in the fiat ribbons is not affected 
by the curvature. The localized nature of the edge states 
remains upon bending, only in some cases the localiza- 
tion lenght increases. In order to illustrate the differences 
between flat and curved geometries the electronic struc- 
ture of the n — 10 zigzag ribbon, both flat and bent, are 
depicted in Figure [8j 

In the absence of SOI, flat and bent ZGNRs only differ 
by a small variation of the subband dispersions. In both 
geometries edge states, n- and cr-orbital derived, appear 
at similar energies although the orbital composition of 
the 7r states is slightly changed with a small admission of 
a orbitals in the curved ribbon. The zero-energy disper- 
sionless states, in the bent geometry, present a very weak 
dispersion. Both, the cr-orbital admission and the disper- 
sion, increase as k increases, and merge with the bulk- 
bands at a k smaller than in the flat geometry. These 
differences are clear in Figure [51 where crossing of the 
edge state occurs in the flat n = 10 ZGNR band structure 
while, due to the hybridization, anticrossing is observable 
in the curved case. SOI effects are similar in both geome- 
tries, although a slight increase of the splittings appears 
in the curved ribbon. 

In order to compare with the n = 10 ribbons, Fig- 
ure [5] represents the band structure of the armchair (5,5) 
carbon nanotube, both without and with SOI. The (5,5) 
CNT, could be thought as formed by joining^ the two 
borders of the curved ZGNR n = 10. Due to the break- 
ing of the rotational symmetry, the number of bands is 



greater in the ribbon than in the tube. While the point 
group of the ZRGN, both flat and bended, is the C3 group 
of graphene, the (5,5) tube belongs to the D 5 d group. 
The corresponding irreducible representations are ten ID 
for the C3 and two ID and four 2D for the D^ 2 ■ Con- 
sequently, considering spin all the subbands are twofold 
degenerated in the ZGNR and two are twofold and four 
fourfold degenerated in the CNT, out of the ten bands 
appearing in the energy interval shown in Figures [8J and 
[5] (see for example Re££). However, in both ribbons and 
NT, the lower bands crossing at the Fermi level are only 
spin degenerate. 

The different degeneracy accounts for the different con- 
sequences that curvature has in the ribbons and in the 
CNTs in the presence of SOI. As it is well known, SOI 
removes all degeneracies compatible with time reversal 
symmetry and parity. Since ZGNR states are only spin- 
degenerated, large splittings comparable to those of the 
fourfold degenerate bands of the NT are not observable. 
Note that the band structure shown in Figure [5] corre- 
sponds to the major possible bended ZGNR, which forms 
an open cylinder with a curvature identical to that of the 
CNT. In fact, curvature effects are also small in ZGNR 
in the presence of the SOI. They reduce to a slight in- 
crease of the energy splitting, see for example the second 
subband. Therefore, enhancement of the SOI strength, 
analogous to that shown for CNT, does not occur in rib- 
bons. 

The effects of curvature on the bearded ribbons are 
very similar, inducing changes in the orbital compositions 
of the edge states as well as a very weak dispersion in the 
flat bands. At T the cr-orbital derived edge zero-energy 
states are shifted down in energy by curvature effects, 
but the shift in energy decreases with the ribbon width. 
Also, there is an small increase of the energy splittings 
in the curved geometry when the SOI is included. 

On the other hand, the main effect of curvature on the 
armchair ribbons is to induce a small gap at T on the 
metallic ribbons between the two linear states crossing 
at the Fermi energy. As well, in both metallic and semi- 
conducting AGNRs, the SOI induced energy splittings 
of the a derived edge states are slightly greater on the 
curved than in the flat geometry. 

Curvature also results in a variation of the spin orien- 
ation axis and in an increase of the orbital angular mo- 
mentum. The expectation value of < S > for the edge 
states is no longer perpendicular to the graphene sheet, 
but always an in-plane contributation turns the spin axis 
to form an angle with the ribbon plane. Further, < L > 
increases by more than two orders of magnitude with 
respect to that of the corresponding flat ribbon. Never- 
theless, for all the ribbons studied it is still much smaller 
than the orbital moment attributed to CNTs due to their 
cylindrical geometry. 

In general, the major effect of curvature is on the edge 
states, in the absence of SOI and for flat ribbons with 
widths greater than 50 chains, the wave functions of the 
degenerate edge states are fully localized at either end of 
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the ribbon, thus in the presence of SOI they become spin- 
filtered states. However, for curved ribbons the states lo- 
calized at the left and right ends interact and accumulates 
at both edges, preventing the formation spin conducting 
chanels. 



IV. DISCUSSION AND CONCLUSIONS 

We have studied the spin-orbit interaction effects on 
the electronic structure of graphene nanoribbons tak- 
ing into account the actual discreteness of the lattice. 
The energy subbands associated to the graphene intrinsic 
electronic bands remain doubly degenerate even when the 
SOI term is included in the Hamiltonian. The atomic lat- 
tice structure satisfies spatial inversion symmetry and, in 
the presence of time reversal symmetry, the intrinsic spin- 
orbit coupling does not break the spin degeneracy^ii. 

The result of the spin-orbit coupling is much more pro- 



nounced in the edge states and particularly in those with 
a large contribution of a orbitals, as it happens in CNTs 
and bilayer graphen e 10 ' 45 . In fact, besides recovering pre- 
vious results on the edge states originated from 7r or- 
bitals, a better understanding of the interplay between 
the intrinsic spin-orbit coupling and the lattice geometry 
is reached. It is found that for the flat zero bands that 
appear at =£- < k < tt and — ^ < k < ^ in the zigzag 
and bearded ribbons respectively, the SOI induced split- 
ting is greater in the regions of the BZ where the larger 
hybridization with a orbitals takes place. 

The inclusion of a orbitals in the basis set gives rise to 
edge states missed in one-band calculations. Edge states 
originated from the cr-orbitals, have been found in the 
three types of ribbons investigated, lying below the Fermi 
energy. These states are localized in the extreme atoms 
at different regions of the BZ (around k = and = ±7r). 
As it occurs for the 7r-derived states, SOI lifts its fourfold 
degeneracy except at the k = ±7r and T points protected 



-1 -0,5 0,5 1 -1 -0,5 0,5 1 
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FIG. 9: Band structure of the armchair (5,5) carbon nanotube calculated with A = (left) and A = OAeV (right). 



by time reversal symmetry. A relevant result is that the 
energy splitting is greater, more than and two orders of 
magnitude, for the a states than for the ir derived states. 
This result confirms those found for CNTs^ and for a 
graphene bilayer where, due to the mixing between tt and 
a bands by interlayer hoppings, the spin-orbit coupling 
is about on order of magnitude larger than in a single 
layer—. 

While the energy spin-orbit coupling induced splittings 
are of the order of 10 -3 meV, equivalent to T« 0.001 K 
for the tt derived edge bands, the edge states originated 
from a orbitals present splittings of the order of 10 _1 
meV i. e. T « 4 - 5 K. These results suggest that a 
states would be particularly appropiatcd to observe the 
spin currents associated with the spin-filtered edge states. 

In the present work, Fermi energy is taken at zero en- 
ergy and neutral graphene is considered, but appropriate 
chemical doping or external gates may change this value 
to tune the energy of the a edge states. In general, they 
lie around 1 eV below the Fermi level, although for the 
Klein GNRs the a derived states are only a few meV away 
from Ep. Therefore, they are easily accesible by apply- 
ing an external potential, which confers Klein nanoribons 
the abbility to exhibit the quantum spin Hall effect. 

Since curvature is known to strongly enhance the spin- 
orbit induced effects on carbon nanotubes, we have con- 
sidered both flat and curved nanoribbons and showed the 
different behaviour that open and close boundary condi- 
tions impose. 

The effects produced by the SOI term and by curva- 
ture would in principle be similar since both induce o-tt 
hybridization and lift degeneracies. Nevertheless, curva- 
ture, for the bearded and armchair ribbons, lifts some 
of the degeneracies at the T point, while SOI does not. 
This result can be understood considering the mapping 



of the two-dimensional graphene bands on the axial di- 
rection of the ribbons: while for the zigzag ribbons the 
Dirac points are projected at k = for armchair K 
and K' points are mapped to T. The curvature of the 
lattice affects the bands around the T region of the BZ in 
the armchair ribbons, in contrast, for the zigzag termina- 
tion, the main effects of curvature occur in the vicinity of 
k = ±^ . Analogously in CNTs curvature induces a gap 
in the primary metal n — 3q zigzag tubes while armchair 
CNTs are metallic. 

The interplay between curvature and SOI effects in 
nanoribbons are not as significant as in CNTs^L The 
boundary conditions imposed by the cylindrical shape 
of CNTs confer rotational symmetry which, alongside to 
the bipartite nature of the honeycomb lattice, is in the 
origin of the unusual electronic properties of CNTs. Or- 
bital magnetic moments 10-20 times larger than the Bohr 
magneton have been observed in suspended CNTs^, at- 
tributed to semiclassical electron orbits encircling the 
tube circumference, whose diameter is much larger than 
the radii of the atomic orbitals^. This behavior does not 
occur in the bent ribbons, where no rotational symme- 
try holds. The spin-orbit coupling effect in bent ribbons 
is mostly due to the increased ir-a hybridization. Thus, 
both curvature and SOI contribute to the increase of the 
orbital angular momentum and to deviate the spin orien- 
tation from the perpendicular direction to the graphene 
sheet. Finally, it is found that curvature does not break 
the chiral symmetry of the edge states and its effects on 
the localization reduce to a small increase of the local- 
ization length. 

In summary, we estimate the magnitude of the intrinsic 
spin-orbit coupling to be of ss 3 meV and although small, 
we have demostrated that, under certain geometries, and 
experimental conditions the Quantum Spin Hall state 
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predicted in graphene should be experimentally observed 
in nanoribbons at temperatures of the order of 3-4 K. 
Further, the spin-filtered edge-states are robust and un- 
affected by curvature. Ours results highlight the impor- 
tance of symmetry to understand the spin-orbit coupling 
effects in graphene nanoribbons. 
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